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In this paper we consider the electric multipole moments of fractal distribution of charges. To de- 
scribe fractal distribution, we use the fractional integrals. The fractional integrals are considered as 
approximations of integrals on fractals. In the paper we compute the electric multipole moments for 
homogeneous fractal distribution of charges. 
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1 Introduction 

Integrals and derivatives of fractional order have found many applications in recent studies in science. 
The interest in fractals and fractional analysis has been growing continually in the last few years. 
Fractional derivatives and integrals have numerous applications: kinetic theories j^EHB]; statistical 
mechanics |1J E] ; dynamics in complex media [71 |H1 El ^1 ; electrodynamics ^1 El El E] and 
many others. The new type of problem has increased rapidly in areas in which the fractal features 
of a process or the medium imposes the necessity of using non-traditional tools in smooth physical 
models. In order to use fractional derivatives and fractional integrals for fractal distribution, we 
must use some continuous medium model [71IS1. We propose to describe the fractal distribution by a 
fractional continuous medium 7^ ^ , where all characteristics and fields are defined everywhere in the 
volume but they follow some generalized equations which are derived by using fractional integrals. In 
many problems the real fractal structure of a medium can be disregarded and the fractal distribution 
can be replaced by some fractional continuous mathematical model. By smoothing of microscopic 
characteristics over the physically infinitesimal volume, we transform the initial fractal distribution 
into a fractional continuous model [3 IE] that uses the fractional integrals. The order of fractional 
integral is equal to the fractal dimension of distribution. The fractional integrals allow us to take 
into account the fractality of the media Fractional integrals are considered as approximations 
of integrals on fractals ^16^. It was proved that integrals on net of fractals can be approximated by 
fractional integrals JBl- In Refs. [H |H1, we proved that fractional integrals can be considered as 
integrals over the space with fractional dimension up to a numerical factor. 

In this paper, we consider electric multipole moments of the fractal distribution of charges. Fractal 
distribution is described by the fractional continuous model [71121111 El- In the general case, the fractal 
distribution cannot be considered as a continuous distribution. There are domains that are not filled 
by particles. We suggest [Z] to consider the fractal distribution as a special (fractional) continuous 



1 



distribution. We use the procedure of replacement of the distribution with fractal dimension by some 
continuous model that uses fractional integrals. The suggested procedure can be considered as a 
fractional generalization of the Christensen approach jl7j . 

In Sec. 2, the density of electric charge for fractal distribution is considered. In Sec. 3, we 
consider the electric multipole expansion. In Sec. 4, the examples of electric dipole moment for fractal 
distribution of chages are derived. In Sec. 5, we consider the electric quadrupole moment of fractal 
distribution of charges. In Sec. 6, the examples of electric quadrupole moments of charged fractal 
parallelepiped are computed. In Sec. 7, the examples of electric quadrupole moments of charged 
fractal ellipsoid are computed. Finally, a short conclusion is given in Sec. 8. 



2 Electric Charge of Fractal Distribution 

Let us consider a fractal distribution of charges. For example, we can assume that charged particles 
with a constant density are distributed over the fractal. In this case, the number of particles N{R) 
enclosed in a volume with characteristic size R satisfies the scaling law 

NiR) ~ i?^, (1) 

whereas for a regular n-dimensional Euclidean object we have N{R) ~ i?". 

For charged particles with number density n{r,t), we have that the charge density can be defined 

by 

p{r,t)=qn{r,t), (2) 

where q is the charge of particle (for example, q is an electron charge), and r = xei + ye2 + ze^. The 
total charge of region W is then given by the integral 

Q{W)=[ p{r,t)dVs=qf n{r,t)dVs, (3) 
Jw Jw 

i.e., Q{W) = qN{W), where N{W) is a number of particles in the region W. The fractional general- 
ization of this equation can be written in the following form: 

Q{W)= [ p{r,t)dVD=q[ n{r,t)dVD, (4) 
Jw Jw 

where D is a fractal dimension of the distribution, and dVn is an element of D-dimensional volume 
such that 

dVD = C3{D,r)dVs. (5) 

For the Riesz definition of the fractional integral 18 , the function C3{D,r) is defined by the 
relation 

The initial points of the fractional integral ^HI are set to zero. The numerical factor in Eq. © 
has this form in order to derive usual integral in the limit D (3 — 0). Note that the usual nu- 
merical factor -f3^{D) = T{l/2)/{2^Tr^/^T{D/2)), which is used in Ref. [IB] leads to 73^^(3 - 0) = 
r(l/2)/(237r3/2r(3/2)) = l/(47r3/2) in the limit D ^ {3 - 0). 
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For the Riemann-Liouville fractional integral ^Sl? the function C^{D^r) is defined by 



Here we use Cartesian's coordinates x, y, and z. In order to have the usual dimensions of the physical 
values, we can use vector r, and coordinates X, i/, z as dimensionless values. 

Note that the interpretation of fractional integration is connected with noninteger dimension |1J 
[21 . This interpretation follows from the well-known formulas for dimensional regularizations. The 
fractional integral can be considered as an integral in the noninteger dimension space up to the 
numerical factor T{D/2)/{2Tr^/'^T{D)). 

If we consider the ball region W = {r : \r\ < R}, Riesz fractional integral and spherically 
symmetric distribution of charged particles {p{r,t) = p{r)), then we have 



^^^^ = r(D/2) io P^''^'' 

For the homogeneous {p{r) = po) fractal distribution, we get 

«<-)^^^«°-"- 

If D = 3, we have Q{W) = {4tt/3)poR^. 

If we consider the Riemann-Liouville fractional integral for the ball region W, and spherically 
symmetric distribution (p{r,t) = p{r)), then we have 

Q(^W) = ^J^^'^Zi^}^^ , r piry'^-'dr. (10) 



r^{D/3)r{D/2) JO 
For the homogeneous = po) fractal distribution, we get 

' DT\D/3)T{D/2) ■ ^ ' 

If D = 3, we get the usual expression Q{W) = {^tt / 2>) p^R^ . 

Fractal distribution of charged particles is called a homogeneous fractal distribution if the power 
law Q{R) ~ R^ does not depend on the translation of the region. The homogeneity property of 
the distribution can be formulated in the following form: For all regions W and W' such that the 
volumes are equal = y(l^'), we have that the number of particles in these regions are equal 

NiyV) = N{W'). Note that the wide class of the fractal media satisfies the homogeneous property. In 
Refs. (3 El) the continuous medium model for the fractal distribution was suggested. 



3 Electric Multipole Expansion for Fractal Distribution of Charges 

It is known that a multipole expansion is a series expansion of the effect produced by a given system in 
terms of an expansion parameter which becomes small as the distance away from the system increases. 
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Therefore, the leading one of the terms in a multipole expansion are generahy the strongest. The first- 
order behavior of the system at large distances can therefore be obtained from the first terms of this 
series, which is generally much easier to compute than the general solution. Multipole expansions 
are most commonly used in problems involving the electric and magnetic fields of charge and current 
distributions, and the propagation of electromagnetic waves. 

To compute one particular case of a multipole expansion, let R = X^ek be the vector from a fixed 
reference point to the observation point; r = Xk^k be the vector from the reference point to a point 
in the distribution; and s = R — r be the vector from a point in the distribution to the observation 
point. The law of cosines then yields 

s"^ = r^ + R'^ -2rR cose (12) 
where s = |s|, r = |r|, R= |R|, and 6 is the polar angle, defined such that 

cos^ = (r,R)/(ri?). (13) 

Using Eq. (|T^ . we get 

s = R]j l-2^cose + ^. (14) 
Now define e = r/R, and ^ = cos 9, then 

But the right hand side of Eq. ()15|) is the generating function for Legendre polynomials -Pn(0 by the 
following relation: 

— 1/2 ^ 

(l-2ee + e2)"' =^e-P„(0, (16) 

n=0 

SO, we have the equation 

n=0 

Any physical potential that obeys a (1/s) law can therefore be expressed as a multipole expansion 

-,00-, 

V=4;i^i:7p«t'-"^»(-%M<'Vz,. (18) 

" n=0 

The n = term of this expansion, called the monopole term, can be pulled out by noting that 
Po{x) = 1, so 



11/" 1 °° 1 /■ 



(18) 



The nth term 

is usually named according to the following: n - multipole, - monopole, 1 - dipole, 2 - quadrupole. 
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4 Electric Dipole Moment of Fractal Distribution of Charges 

An electric multipole expansion is a determination of the voltage V due to a collection of charges 
obtained by performing a multipole expansion. This corresponds to a series expansion of the charge 
density p{r) in terms of its moments, normalized by the distance to a point R far from the charge 
distribution. In MKS, the electric multipole expansion is given by Eq. ((TH|l: 

V=^,^J^ir"P„(cosmrWo, (21) 

where Pn{cos9) is a Legendre polynomial. 

The first term arises from PoiS,) = 1, while all further terms vanish as a result of PniC) being a 
polynomial in ^ for n > 1, giving Pn{0) = for all n > 1. 

Set up the coordinate system so that 9 measures the angle from the charge-charge line with the 
midpoint of this line being the origin. Then the n = 1 term is given by 

^1 = I rPi{cose)p{r)dVD = 

47reo P- Jw 



D 



47reoi2^ 

For a continuous charge distribution, the electric dipole moment is given by 



pW = / rp{r)dV3, (23) 
Jw 

where r points from positive to negative. Defining the dipole moment for the fractal distribution by 
the equation 

p(^) = / rp{r)dVD, (24) 
Jw 

then gives 

^ 1 (p(^),R) ^ 1 p(^)cosa 
' 47reo 47r£o R2 ' ^^^^ 

where we use 



cosa = (p(^),R)/(p(^)i?), p(^) = V(prO' + (4 T + (Pzf- (26) 

Let us consider the example of electric dipole moment for the homogeneous (/o(r) = po) fractal 
distribution of electric charges in the parallelepiped region W: 



0<x <A, 0<y<B, 0<z<C. (27) 
In the case of Riemann-Liouville fractional integral, we have pi^^ in the form 

T-^[a) Jq Jo Jq T («) a^(a + l) 
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where a = D/2>. The electric charge of parallelepiped region (|27() is defined by 



Therefore, we have the dipole moments for fractal distribution in parallelepiped in the form 

P?^ = -^^Q{W)A, (30) 
a + 1 



By analogy with these equation, we can derive 

P^^^ = ^Q{W)B, pP) = ^Q(t^)a (31) 



Using a/{a + 1) = D/{D + 3), we get 



p'''' = ^/'\ (32) 



where p^^^ = |p^^^| are the dipole moment for the usual 3-dimensional homogeneous distribution. For 
example, the relation 2 < < 3 leads us to the following inequality 

0.8 < p(^Vp^^^ < 1- (33) 



5 Electric Quadrupole Moment of Fractal Distribution of Charges 

There are also higher-order terms in the multipole expansion that become smaller as R becomes large. 
The electric quadrupole term in MKS is given by 

^2 = -r—^ I r^P2{cose)p{r)dVD = 

The electric quadrupole is the third term in an electric multipole expansion, and can be defined 
in MKS by 

where Eq is the permittivity of free space, R is the distance from the fractal distribution of charges, and 
Qki is the electric quadrupole moment, which is a tensor. Note that are Cartesian's coordinates 
of the vector R, and Xk are coordinates of the vector r. 

The electric quadrupole moment is defined by the equation 

Qki= [SxkXi - r'^6M]p{r)dVD, (36) 
Jw 
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where = x,y, or z. From this definition, it fohows that 



3 

Qki = Qik, and Qkk = 0. (37) 
fc=i 



Therefore, we have Qzz = —Qxx — Qyy In order to compute the values 



Qi^,> = / [2x' -y'- z']p{r)dVD, (38) 
Jw 

Qiy^ = Jj-^^ + V - z^W)dVD, (39) 

Qif ^ = / [-^' -y' + 2z^]pir)dVD, (40) 
Jw 

we consider the following expression 

Q{a,l3n)=f [ax^ + Py^ + 7z^]p{r)dVD. (41) 
Jw 

Using Eq. (|1T|). we have 

Q(^)=Q(2,-1,-1), Q(^)=Q(-1,2,-1), Q(f) = Q(-l,-l,2). (42) 

The example of electric quadrupole moment for the parallelepiped and ellipsoid regions are considered 
in Sees. 6 and 7. 

6 Quadrupole Moment of Charged Fractal Parallelepiped 

Let us consider the example of electric quadrupole moment for the homogeneous (/o(r) = po) fractal 
distribution of electric charges in the parallelepiped region ()27() . If we consider the region W in the 
form (P7|). then we get 

[a + 2ja^i ■^[a) 

where we use the Riemann-Liouville fractional integral '18', and the function 0^(0, r) in the form 

C3(Ar) = ^|i— , a = D/3. (44) 
The electric charge of the region W is 

If I? = 3, we have Q{W) = poABC. Using equations and we get the following equation 

Q(q,/3,7) = ^Q(T^)M2 + + ^C2]. (46) 
a + 2 
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If D = 3, then we have a/{a + 2) = 1/3. As the result, we have electric quadrupole moments Q^^'' of 
fractal distribution in the region W: 

(D) _ 3£> (3) 
^kk ~ D j^Q ^kk J ) 

where Qj^j! are moments for the usual homogeneous distribution {D = 3). By analogy with these 
equations, we can derive Q^.^ ^ for the case k ^ I. These electric quadrupole moments are 

Q'S^ = J^TW ' (MO- (48) 
Using inequality 2 < D < 3, we get the relations for diagonal elements 

0.75 < qI^VqS < 1, (49) 



and nondiagonal elements 
where k ^ I. 



0.64 < Qif M? < 1- (50) 



7 Quadrupole Moment of Charged Fractal Ellipsoid 

Let us consider the example of electric quadrupole moment for the homogeneous (/o(r) = po) fractal 
distribution in the ellipsoid region W: 

^ + ^ + ^<i. (51) 

If we consider the region W in the form H51|). then we get expression (|41() in the form 

= r^^'llfr^)"^ M^Zi(a) + (3B'Z2{a) + jC^Zsia)], (52) 
(3o + 2)T-^{a) 

where a = D/3, and Zi{a), i = 1,2,3 are defined by 

Zi{a) = S{a + l,a- l)S{a - 1, 2a + 1), (53) 

^2(0) = 5(a - 1, a + l)S{a - 1, 2a + 1), (54) 
Z^ia) = S{a-l,a- l)S'(a + 1,2a - 1). (55) 
Here we use the following function 

^ N ■ N r(n/2 + l/2)r(m/2 + 1/2) 

S(n,m) = / dx cos"'(x) sm"'(x) = — ; — / (56) 

^ ^ Jo 2r(n/2 + m/2 + 1) ^ ^ 



Note that Zi{a) = Z2{a) = .^3 (a). Using these equations, we get the following relation: 

Po(^^C)°r^(a/2)r(a/2 + ] 
(3a + 2) r3(a)r(3a/2 + 1 
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Using r(/3 + 1) = /?r(/?), we have 



If = 3, we obtain 



Q(a,A7) = ^^^^#^M^ + /3i?^ + 7C^^]. (59) 
15 



The total charge of the elhpsoid region W is defined by 

Q(W) = p„ = P.(.4«C)° 3„r3(!)r(3i/2) - <^°> 

If L> = 3, we have the total charge Q(iy) = {4Tr/3)poABC . Here we use r(l/2) = 

Using Eqs. and (jHOj), we can derive the electric quadrupole moments for fractal ellipsoid. 
As the result, we have 

Q{c^,P,l) = in^QiW)[aA' + (3B^ + jC% (61) 

dJJ + D 

If I? = 3, then we have the well-known relation: 

Q{a,(3,-f) = {l/5)Q{W)[aA^ + (3B^ + jC% (62) 

If 2 < Z) < 3, then we have 

I < Q'^/Q^l < 1- (63) 

The nondiagonal elements of electric quadrupole moment are defined by the following equations 

Qxy = 3po / xy dVo, Qxz = 3/90 / xz dVn, Qyz = 3po / yz dVo- (64) 
Jw Jw Jw 

Using these equations, we can derive the nondiagonal elements in the form: 



3a + 2 


r3(a)r(3a/2 + l) 


Qpo{ABCY 


r(a/2)r2(a/2 + l/2) 


3a + 2 


r3(a)r(3a/2 + 1) 


QPq{ABCT 


r(a/2)r2(a/2 + l/2) 


3a + 2 


r3(a)r(3a/2 + l) 



Using r(/3 + 1) = j3T{f3) and Eq. (|6U|) . we get the following equations 
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As the result, we have 

where k ^l. Here we use T{l/2) = -y/vr 

0.6972 < Q'^^/Qti' < 1. (72) 

8 Conclusion 

In this paper, we use fractional continuous model for fractal distribution of electric charges. The 
fractional continuous models for fractal distribution of particles can have a wide application. This is 
due in part to the relatively small number of parameters that define a fractal distribution of great 
complexity and rich structure. In many cases, the real fractal structure of matter can be disregarded 
and the distribution of particles can be replaced by some fractional continuous model [3 13- In order to 
describe the distribution with noninteger dimension, we must use the fractional calculus. Smoothing 
of the microscopic characteristics over the physically infinitesimal volume transforms the initial fractal 
distribution into fractional continuous model that uses the fractional integrals. The order of fractional 
integral is equal to the fractal dimension of the distribution. The fractional continuous model for 
the fractal distribution allows us to describe dynamics of a wide class fractal media [HI El 01 • One 
of the dynamical equation of physics is a Liouville equation. Note that the Liouville equation is a 
cornerstone of the statistical mechanics. The fractional generalization of the Liouville equation was 
suggested in Refs. jHEl- The fractional generalization of the Liouville equation allows us to derive 
the fractional generalization of the Bogoliubov equations Using fractional analog of the Liouville 
equation and Bogoliubov equations OiOl j we can derive the description of a fractal distribution 
as a fractional system. Fractional systems can be considered as a special case of non-Hamiltonian 
systems jH j^. Note that non-Hamiltonian systems can have stationary states of the Hamiltonian 
systems [H HOI iH (111113121 . 
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